Non-equilibrium structures and slow dynamics in a two dimensional spin system with 
competitive long range and short range interactions 
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We introduce a lattice spin model that mimics a system of interacting particles through a short 
range repulsive potential and a long range attractive power law decaying potential. We perform 
a detailed analysis of the general equilibrium phase diagram of the model at finite temperature, 
showing that the only possible equilibrium phases are the ferromagnetic and the antiferromagnetic 
ones. We then study the non equilibrium behavior of the model after a quench to subcritical tem- 
peratures, in the antiferromagnetic region of the phase diagram region, where the pair interaction 
potential behaves in the same qualitative way as in a Lennard- Jones gas. We find that, even in 
the absence of quenched disorder or geometric frustration, the competition between interactions 
gives rise to non-equilibrium disordered structures at low enough temperatures that strongly slow 
down the relaxation of the system. This non-equilibrium state presents several features character- 
istic of glassy systems, such as subaging, non trivial Fluctuation Dissipation relations and possible 
logarithmic growth of free energy barriers to coarsening. 
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I. INTRODUCTION 

The nature of glassy magnetic states in the absence 
of quenched disorder has been object of a great deal 
of work!, both experimental and theoretical. Simple 
experimental realizations of nondisordered systems in 
which glassy phases have been found are antiferromag- 
nets (AFM's) with kagome geometries^^. In particu- 
lar the presence of slow dynamics and aging effects in 
kagome antiferromagnets is well established—. Anyway, 
these are non-disordered geometrically frustrated sys- 
tems, so the understanding of the mechanisms present 
on the dynamics of these systems has to deal with the 
effects of the involved geometry of the kagome lattice. 
On the other hand, glassy behavior in structural glasses 
appears dynamically, without any kind of imposed dis- 
order or geometrical frustration. A prototype model for 
structural glasses is the Lennard-Jones binary mixture^. 
One may wonder whether glassy behavior can appear in 
lattice spin systems sharing some of the basic features 
of the Lennard-Jones model, such as the competition be- 
tween short range repulsive interactions (i.e., hard core) 
and long range attractive interactions. A simple model 
with those properties is the Ising model with competitive 
interactions on the square lattice. 

Consider the general lattice Hamiltonian 

H-J^w + J^y- (1) 

(id) (id) 13 

where a = ±1. The first sum runs over all pairs of nearest 
neighbor spins on a square lattice and the second one over 
all distinct pairs of spins of the lattice, is the distance, 
measured in crystal units between sites i and j. For 



Ji < and J2 > this Hamiltonian describes an ultra- 
thin magnetic film with perpendicular anisotropy in the 
monolayer limh>£ and it has been the subject of several 
theoretical studies (see Refsl^H and references therein). 

In this work we consider the case J2 = — 1 (long range 
ferromagnetic interactions) J = J\j \ J%\ > (short range 
antiferromagnetic interactions), so (UJ) reduces to the di- 
mensionless Hamiltonian 

(id) (id) 13 

For J > 1 this Hamiltonian mimics a system of parti- 
cles interacting through a short range repulsive potential 
and a long range power law decaying attractive poten- 
tial, qualitatively similar to the Lennard-Jones pair in- 
teractions potential (see an example in FigQ]). Although 
the exponent 3 in the power law interacting potential is 
arbitrary in this case, it presents two advantages. First, 
in a two dimensional system it is large enough to ensure 
the existence of the thermodynamical limit and, second, 
it allows us to use all the previous knowledge of the much 
more studied related system J < and J\ > (ultrathin 
magnetic films model). As we will show, it displays com- 
plex low temperature dynamical behavior, even when its 
equilibrium properties are simpler than those observed 
in the ultrathin magnetic films case. In order to corre- 
late equilibrium and non-equilibrium properties we start 
our analysis by investigating the finite temperature ther- 
modynamical behavior of the model. Since, to the best 
of our knowledge, this model has not been previously 
studied in the literature and also for completeness we 
perform in section [H] a detailed analysis of the complete 
equilibrium phase diagram using Monte Carlo (MC) sim- 
ulations. In section UlTl we analyze the low temperature 



2 



0.5 
0.4 
0.3 
*C 0.2 

-©- 0.1 

0.0 
-0.1 
-0.2 




4 6 

r 



10 



FIG. 1: Pair interaction potential as a function of distance 
for J = 1.5. 



relaxation properties of the model, by computing differ- 
ent quantities like the average linear size of domains, en- 
ergy, two times correlation functions and Fluctuation- 
Dissipation relations. In section [IV] we discuss our re- 
sults, comparing them with previous reported results of 
slow dynamics in non-disordered systems, in particular, 
the Lennard- Jones gas. 



II. EQUILIBRIUM PHASE DIAGRAM 

As a first step we analyze the zero temperature prop- 
erties of the model. In the absence of the long range term 
the model reduces to the Ising antiferromagnetic model 
and the ground state of the system is a Neel antiferromag- 
netic state. For J = it is clear that the ground state 
is ferromagnetic. To obtain the ground state between 
these two limits we evaluated the energy per site for dif- 
ferent spin configurations, namely ferromagnetic, anti- 
ferromagnetic, stripes of different widths (1,2,... rows 
of spins) and chequered domains of different sizes. The 
energies per spin of the ferromagnetic and antiferromag- 
netic states are given by Ef = 2 J — a and E a f = —2J + b 
respectively, with±£ a = 4.5168 and b = 1.3230. The 
energy per spin of a state composed by ferromagnetic 
stripes of width h is given by E s {h) = 2 (1 — 1/h) J + Sh- 
The values of Sh were calculated numerically in Refflpl: 
for instance Si = 0.4677, S2 = —0.7908, etc.. In Figure 
[2] we compare the energy per site for different configura- 
tions. The figure shows that the only stable states are 
the ferromagnetic and the antiferromagnetic ones for any 
value of J. By equating Ef = E a f we obtain for the tran- 
sition point between the ferromagnetic and the antiferro- 
magnetic states the value J t = 1.4599: the ground state 
is ferromagnetic below this value and antiferromagnetic 
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FIG. 2: (Color on-line) Energy per site vs J for different spin 
configuration; hi corresponds to an state composed by stripes 
with width one, h2 to stripes with width two and so on. The 
chequered states give similar values for the energy per site 
(not shown here) to the stripes states with the same width. 



above of it. We also checked different chequered anti- 
ferromagnetic states^, verifying that they have higher 
energies than either the ferromagnetic or the antiferro- 
magnetic states for any value of J. Monte Carlo (MC) 
simulations at low temperatures confirm that these are 
the only low temperature stable phases. 

We next consider the equilibrium finite temperature 
properties of the model, by using Metropolis Monte Carlo 
algorithm in finite square lattices with N = L x L sites 
and periodic boundary conditions. To handle the contri- 
bution of the long range terms in the periodic boundary 
conditions we used the Ewald sums technique^. In order 
to speed up the simulations, the codes were implemented 
by keeping track of all the local fields. In this way, local 
fields update (an operation which is of O(N)) is per- 
formed only when a spin flip is accepted. This imple- 
mentation is very effective when relaxation is very slow 
(in general at low temperatures) and therefore the accep- 
tance rate is small, while it does not change the computa- 
tional time when the acceptance rate is high (usually at 
intermediate or high temperatures). This will be particu- 
larly important when considering non equilibrium effects 
in section IIII1 allowing us to treat large system sizes at 
very low temperatures. 

To characterize the critical properties of the model we 
calculate different thermodynamical quantities as a func- 
tion of the temperature, for different values of J and 
L, namely, the magnetization per spin m, the staggered 
magnetization per spin m s , the associated susceptibilities 



X (T) = ^«m 2 } - (m) 2 ) 



(3) 
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Xs(T) = ^((m 2 s ) 



(m s ) 2 ), 



the specific heat 



C(T) = j^((H 2 } (if) 2 ) 



(4) 



(5) 



and the fourth order cumulant 



3(H 2 ) 2 
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where (. . .} stands for an average over the thermal noise. 
All these quantities are calculated starting from an ini- 
tially equilibrated high temperature configuration and 
slowly decreasing the temperature. For every temper- 
ature the initial spin configuration is taken as the final 
configuration of the previous temperature; we let the sys- 
tem to equilibrate Mi Monte Carlo Steps (one MCS is 
defined as a complete cycle of N spin update trials) and 
average out the results of M 2 MCS, typical values of Mi 
and M2 being around 10 5 and 5 x 10 5 respectively. 

Figures Q and Q show the typical results for the dif- 
ferent thermodynamical quantities in the antiferromag- 
netic region of the phase diagram, namely, for J > J t . 
Fig.(@jD) shows that the fourth order cumulant exhibits a 
vanishing minimum, consistent with a second order phase 
transition. Fig. (J3Jd) shows that the staggered suscepti- 
bility exhibits a size dependent maximum which scales 
as L 1 ^, with j/v — 1.8 ± 0.1, consistent with the ex- 
act value 7/V = 1.75 of the two dimensional short range 
Ising model. Moreover, as J increases j/v approaches 
systematically the value 1.75 (for instance, for J = 3 
we found y/v = 1.74 ± 0.05; see Fig. ((TJ) . We see from 
Fig.(j4^i) that the specific heat exhibits a size dependent 
maximum which scales as L a ' v , with a/v = 0.23 ± 0.05; 
similar values were found for other values of J > Jt (see 
FigtZj) . Although small, those values are larger than ex- 
pected for a phase transition in the universality class of 
the two dimensional Ising model (a = 0). However, since 
those values are also observed for large values of J, we 
believe that this is a finite size effect. Hence, we con- 
clude that the whole line between the paramagnetic and 
the antiferromagnetic phases (J > J t ) belongs to the uni- 
versality class of the short range two dimensional Ising 
model. 

The critical properties for J < Jt are a bit more com- 
plex. For J < 1.3 the order parameter (magnetization), 
the susceptibility, the specific heat and the fourth or- 
der cumulant present qualitatively the same behavior as 
those quantities in the J > Jt case, but with a different 
set of critical exponents (see Fig[7j). We found y/v w 1.1 
and a/v « 0.14., which are close to the renormalization 
group estimates for the J = caseii: v = 1, a = 
and 7 = 1; the small difference between those values 
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FIG. 3: (Color on-line) Order parameter (absolute value of the 
staggered magnetization) (a) and associated staggered suscep- 
tibility (b) as a function of the temperature for J = 1.6 and 
for different system sizes; the inset in (b) shows the finite size 
scaling of the maximum of \s- 



and ours can be attributed to finite size effects, which 
are very strong when the long range ferromagnetic in- 
teractions dominate. Hence, we conclude that the whole 
line for < J < 1.3 belongs to the universality class of 
the two dimensional 1/r 3 ferromagnetic Ising model. For 
1.3 < J < Jt we find a clear evidence that the ferro-para 
transition is a first order one. The typical behavior of 
the thermodynamical quantities in this case is illustrated 
in Figs. [5] and 03 We see that the fourth order cumulant 
presents a clear converging minimum as the system sizes 
increases, as expected in a first order transition^. The 
finite size scaling of susceptibility is also consistent with 
the L 2 behavior expected for a first order transition in 
a two dimensional system^. The specific heat exponent 
for J = 1.4 is a/v = 1.1 ± 0.2. This value is certainly far 
from 2 (the expected value in a first order transition) , but 
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FIG. 4: (Color on-line) Moments of the energy as a function 
of the temperature for J = 1.6 and for different system sizes, 
(a) Specific heat C; the inset shows the finite size scaling of 
the maximum of C. (b) Fourth order cumulant. 
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FIG. 5: (Color on-line) Order parameter (absolute value of 
the magnetization) (a) and associated susceptibility (b) as 
a function of the temperature for J = 1.4 and for different 
system sizes; the inset in (b) shows the finite size scaling of 
the maximum of %■ 



it is larger than the critical exponent of any continuous 
transition. Besides finite size effects, such large difference 
is probably also associated to the presence of a tricritical 
point somewhere between J = 1.3 and J = 1.4. This as- 
sumption is consistent with the fact that ajv approaches 
the expected value a/v — 2 as J increases approaching 
J = Jt (we obtained ajv — 1.8 ± 0.1 for J = 1.43; see 
FigEJ. 

We summarize the obtained results for the critical ex- 
ponents in FigE] and the overall phase diagram in FigEl 



III. NON-EQUILIBRIUM PROPERTIES 

This section deals with the far-from equilibrium prop- 
erties of the system at low temperatures, i.e., its relax- 
ation dynamics after a sudden quench from T = oo to a 



temperature T < T c . 



A. Non equilibrium domain structures: energy 
relaxation and characteristic domain length 

First we analyze the time evolution of the energy, 
with the time measured in MCS. We consider both the 
instantaneous energy per spin E/N (i.e., the energy 
along single MC runs) and the mean excess of energy 
Se(t) = [H]/N — u(T), where [. ..] stands for average 
over different MC runs (i.e., over different realizations of 
the thermal noise). u{T) is the equilibrium energy per 
spin at temperature T; u(T) is obtained by equilibrating 
first the system during 10 4 MCS starting from the ground 
state configuration and then averaging over a single MC 
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FIG. 6: (Color on-line) Moments of the energy as a function 
of the temperature for J = 1.4 and for different system sizes, 
(a) Specific heat C; the inset shows the finite size scaling of 
the maximum of C. (b) Fourth order cumulant. 



run during 10 5 MCS. 

To check out our results we first calculate the evolu- 
tion of Se(t) in the simple case J < J t for different quench 
temperatures. The typical behavior is shown in Figure 
[9l We find that, after a short transient period and be- 
fore the system completely relaxes, the excess of energy 
behaves as Se(t) ~ t~ x / 2 independently of T. Since it 
is expected that Se(t) oc l/l(t), where l(t) is the char- 
acteristic length scale of the domains, this behavior is 
consistent with a normal coarsening process of a system 
with non-conserved order parameter—, where l(t) ~ t 1 ! 2 . 

Next, we consider the relaxation in the antiferromag- 
netic region J > Jt for different quench temperatures. At 
low enough temperatrues the relaxation of the system 
clearly departs from that expected in a normal coars- 
ening process. The typical behavior of the instanta- 
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FIG. 7: Critical exponents obtained from finite size scaling as 
a function of J. (a) Susceptibility exponent 'y/u; (b) Specific 
heat exponent a/ v. The dashed lines indicate the reference 
values 1, 1.75 and 2 in (a) and and 2 in (b). 



neous energy is shown in Fig[l0l together with typical 
domain configurations along single MC runs for J = 2 
and T — 0.04. In that figure the domains correspond to 
regions of antiferromagnetic ordering, namely, black and 
white colors codify regions with local staggered magneti- 
zation m s w 1 and m s w — 1 respectively. 

Different relaxation regimes can be identified. After 
a short time quick relaxation process < t < tq w 
20 MCS, in which local antiferromagnetic order is set, 
the system always gets stuck in a complex non equilib- 
rium disordered state composed mainly by a few inter- 
mingled macroscopic antiferromagnetic domains; its typ- 
ical shape is illustrated for a larger system size in FigfTTl 
This state presents a sort of labyrinth structure, in the 
sense that there is always at least one macroscopic con- 
nected domain, i.e., in such domain any pair of points 
can be connected by a continuous path without cross- 
ing a domain wall (see for example the black domain in 
Fig fTTl) . Up to certain characteristic time t\ the system 
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FIG. 8: Phase diagram T vs. J. The critical temperatures 
were estimated from the maxima of the specific heat. Filled 
circles and open hexagons correspond to second and first order 
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FIG. 10: Instantaneous energy per spin as a function of time 
for single realizations of the stochastic noise for J = 2, L = 48 
and T = 0.04. Typical antiferromagnetic domain config- 
urations are shown along the evolutions, where black and 
white colors codify regions with local staggered magnetiza- 
tion m s ?s 1 and m s fa — 1 respectively. The dashed lines 
correspond to the equilibrium energy at this temperature, (a) 
After living the glassy regime the system equilibrates, (b) Af- 
ter living the glassy regime the system gets stuck in a striped 
configuration. 



FIG. 9: (Color on-line) Excess of energy 5e(t) as a function 
of time for J = 1, L = 100 and different quench temperatures 
T < T c . The results were averaged over 2000 MC runs. 

slowly relaxes by eliminating small domains and fluctua- 
tions located in the large domain borders, in such a way 
that the local curvature of the domain walls is reduced 
(FigHDf . Along this process the area of the main domains 
remains almost constant. In this sense, such process is 
reminiscent of a spinodal decomposition. We will call this 
the glassy regime. For time scales longer than a certain 
characteristic time t\ both domains finally disentangle 
and relaxation is dominated by the competition between 
only two large domains, separated by rather smooth do- 
main walls. Fig[10] illustrates the two possible outcomes 
of this process: either the system relaxes directly to its 
equilibrium state (FigJIDk) or it gets stuck in an ordered 
configuration composed of stripe shaped antiferromag- 
netic domains with almost flat domain walls (FigfTUb). 




FIG. 11: Typical antiferromagnetic domain configuration for 
J = 2, T = 0.04, L = 256 and t = 100; MCS. Black and 
white follows the same convention as in FigllOl 
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FIG. 12: (Color on-line) Excess of energy per spin (see text 
for details) as a function of time for J = 2, T = 0.06 and 
different system sizes. Every curve was obtained by averaging 
over 400 runs. The dashed and full lines correspond to a power 
law and exponential fittings respectively. 



FIG. 13: (Color on-line) . Excess of energy per spin as a 
function of time for J = 2, L = 48 and different quench 
temperatures T < T c (decreasing from left to right). Every 
curve was obtained by averaging over 400 runs. The inset 
shows an Arrhenius plot of the crossover time. ri. 



We observe that both outcomes can happen with finite 
probabilities, the former being a bit more probable than 
the latter. The second case covers a large variety of con- 
figurations, including more than two stripes that can be 
oriented parallel to one of the coordinate axes (as in 
FigfTOb) or diagonally oriented (not shown). We will 
call this the ordered regime. Once the system arrives to 
one striped configuration, relaxation proceeds through 
the parallel movement of the domain walls, which per- 
form a sort of random walk until two walls collapse and 
the system either attains the equilibrium state or gets 
stuck in a new striped configuration with a lesser number 
of stripes. The mechanism of movement of the domain 
walls in this case is dent formation, i.e., single isolated 
spin flips along the interface creating an excitation that 
propagates along it, until either it disappears or covers 
the whole line^, which therefore advances in the per- 
pendicular direction. The same kind of non-equilibrium 
structures and relaxation dynamics has been observed 
in two dimensional short range interacting spin models 
at very low temperatures, namely the Ising^ or Potts^ 
models. However, in those cases the movement of the 
dents are dominated by single spin flip barriers, while in 
the present one the associated mechanism is more com- 
plex due to the long range interactions. 

In Fig[12] we illustrate the typical behavior of the ex- 
cess of energy 5e(t) at a fixed temperate for different 
system sizes. In FigH3]we show the excess of energy for 
different temperatures at a fixed system size. The three 
different relaxation regimes can be clearly seen in those 
curves: transient, glassy and ordered. The glassy regime 
appears for temperatures smaller than certain value T g 
(T g w 0.15 for J = 2). In this regime the excess of en- 
ergy exhibits a size-independent pseudo-plateau, where 



it decays very slowly; indeed, the behavior of Se(t) can be 
well fitted by a power law Se(t) ~ with very small 
exponents that decrease with temperature (the exponent 
for J = 2 ranges from oj « 0.03 for T=0.04 up to w « 0.1 
for T = 0.1), suggesting a logarithmic relaxation at very 
low temperatures. This suggests an activated dynam- 
ics with multiple energy barriers (we will return to this 
point later). After this regime, the system relaxes expo- 
nentially into the ordered regime Se(t) ~ e - t l T ^ T ) (see 
Fig fT2|) . The characteristic relaxation time n(T) can be 
estimated by fitting the corresponding part of the relax- 
ation curve, as shown in Fig|T2j The inset of Fig fl3l shows 
an Arrhenius plot of n . The exponential decay, together 
with the clear Arrhenius behavior of n, indicates that 
the crossover between the two regimes is dominated by 
the activation through a single free energy barrier. 

To gain further insight about the nature of the relax- 
ation in the glassy regime, we analyze the scaling prop- 
erties of the characteristic domain length l(t). A sensible 
way to estimate the behavior of that quantity is to define 



it as 



17.18.19 



l(t) 



-u(T) 
Se(t) 



(7) 



In FigfUwe show l(t) for J = 2, L = 48 and differ- 
ent temperatures T < T g . The behavior of the excess of 
energy implies that, for time scales to < t < n, lit) 
increases very slowly from a temperature-independent 
value lo — 1(tq); for time scales t > t\ the charac- 
teristic length departs exponentially from the pseudo 
plateau (see Figfl4k). We estimated Z as the average 
of the curves for different temperatures at r , obtaining 
Iq w 5.78. In FigfT4b we show a double log plot of the 
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FIG. 14: (Color on-line) (a) Characteristic domain length 
(see text for details) as a function of time for J = 2, L = 48 
and different temperatures T < T c (decreasing from left to 
right), (b) Log-log plot of the normalized length (l — lo)/T a vs. 
ln(t) from the same data as in (a) for the lowest temperatures; 
lo — 5.78 is indicated in (a); the value of the exponent a = 2.7 
was chosen to obtain the best data collapse of the curves in 
the glassy regime. The straight line is a reference (power with 
exponent a). 



rescaled quantity (l(t) — lo)/T a vs. ln(t). The expo- 
nent a was chosen to obtain the best data collapse in the 
glassy regime of the data presented in FigfT4k. Actually, 
a good data collapse inside the error bars of the sta- 
tistical fluctuations is obtained for values of a between 
2.65 and 2.75; for values of the exponent outside that 
range the curves clearly do not collapse. Hence, we esti- 
mated a = 2.7 ±0.05. The power law like behavior of the 
rescaled curves in Figfl4b shows that the characteristic 
length behaves as 
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FIG. 15: (Color on-line) Arrhenius plot of the characteristic 
time for shrinking squares for different values of the square 
side l s : from bottom to top l s = 5,6,7,8,9. The lines are a 
guide to the eye. The inset shows the energy barrier (slope 
of the linear fittings in the Arrhenius plot) as a function of 
l s (the symbol size in the inset is larger than the statistical 
error bars). 



for to < t < ri(T) (a log-log plot of l(t) — Iq vs. t shows 
that a power law fit in the entire time interval is clearly 
inferior than in Fig fT4|) . Such behavior is consistent with 
a class 4 system, according to Lai et al classification 2 ^, 
i.e., a system with domain size dependent free energy 
barriers to coarsening^ f{l). In our case this would cor- 
respond to f(l) ~ b(l — lo) 1 ^- The numerical results 
suggest that in the present model such growth would 
stop when some maximum characteristic length l ma x is 
reached at n(T), where the barrier becomes indepen- 
dent of I. After this point the system relaxes exponen- 
tially with a characteristic time t\ oc exp(F/T), where 
F = f{lmax) and therefore l max ~ Iq + {F/b) a . From 
the data of Figfl4b we estimate b « 0.32, while from the 
data of the inset of Fig[13]we estimated F ~ 0.44, giving 
an estimation l m ax ~ 8. 

To check the above interpretation we analyze the char- 
acteristic time t s for shrinking squares, i.e., the time 
needed for a square excitation of linear size l s to com- 
pletely relax. This technique has been proved to be a sen- 
sitive way to check the relaxation dynamics of short range 
models when free energy barriers are involve d 18 ' 19 ' 21 . In 
particular, Shore et aU£ have argued (and shown to be 
valid in particular cases) that the energy barriers to 
shrink square-shaped excitations should be a measure of 
the free energy barriers to coarsening. In our case we 
started with a ground state configuration of size L with 
a square of inverted spins of size l s (L ^ l s ) and peri- 
odic boundary conditions. Although it is not clear to us 
whether the arguments of Shore et aU^ can be straightfor- 
wardly extended to a system with long range interactions 
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or not, one can still expect the barriers to shrink a square 
to provide at least a rough measure of the free energy bar- 
riers to coarsening. In our case, this expectation is based 
on the direct observation of the domain configurations 
during relaxation in the glassy regime. We observe that 
rough domain walls tend to become flat rather fast, and 
that relaxation proceeds mainly at small jumps in the 
energy every time a sharp edge moves. The results for 
the time for shrinking squares support this conjecture. 
In Fig[l5] we show an Arrhenius plot of t s for different 
values of l s and temperatures T <T g . We see that t s ex- 
hibits a clear Arrhenius behavior at all the temperatures 
for l s > 4 (for sizes l s < 4 the squares shrink quickly in 
a few MCS), with associated barriers that grow slowly 
for l s < 8 and saturate for l s > 9 at a value around 0.5, 
close to F = 0.44. Although the limited range of values 
of l s where the barrier shows a dependency on it does 
not allow a more accurate comparison, the consistency 
with the previous interpretation of the behavior of l(t) is 
clear. 

For temperatures larger than T g the glassy regime com- 
pletely disappears and the system decays through a nor- 
mal coarsening process, i.e., Se(t) ~ i -1 / 2 (see FigJT3]) . 
However, for some range of temperatures it still gets 
stuck in some long-lasting antiferromagnetic striped con- 
figuration with high probability, so the corresponding 
plateau in the excess of energy is still observable (for 
J = 2 we observed it for temperatures up to T ~ 1.5). 
Those configurations are highly stable, even at relatively 
high temperatures. The characteristic equilibration time 
T2 , defined as the time after which the system attains the 
equilibrium state with probability one, is very difficult 
to estimate, but it is at least three orders of magnitude 
larger than n for T <T g . 
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FIG. 16: (Color on-line) Two-times autocorrelation function 

C(t w + t,t w ) as a function of (t + t w )/t w for J = 1 (ferro- 
magnetic phase), T=0.2, L = 600 and different values of the 
waiting time t w (increasing from bottom to top). 



Although aging can be detected through different time- 
dependent quantities, a straightforward way to establish 
it in a numerical simulation is to calculate the spin auto- 
correlation function 



1 N 

^ &i{tw + t)(Ti(t w ) 
i=l 



(9) 



B. Time correlation and response functions 

Another way to characterize the out of equilibrium dy- 
namics of complex magnetic systems is through the anal- 
ysis of the two-time autocorrelation function C(t,t'). A 
system that has attained thermodynamical equilibrium 
or meta-equilibrium satisfies time translational invari- 
ance (TTI), i.e., C(t,t') = C{t - t'), al least for certain 
time scales. Far from equilibrium TTI is broken and 
time correlations exhibits a dependency on the history of 
the sample after the quench. This phenomenon is called 
aging and in real systems it can be observed through a 
variety of experiments. A typical example is the zero- 
field-cooling^ 2 - experiment, in which the sample is cooled 
in zero field to a subcritical temperature at time t = 0. 
After a waiting time t w a small constant magnetic field 
is applied and the time evolution of the magnetization is 
recorded. It is then observed that the longer the waiting 
time t w the slower the relaxation and this is the origin 
of the term aging. Moreover, the scaling properties of 
two-times quantities provides information about the un- 
derlying relaxation dynamic a 23 i 24 . 



where t w is the waiting time from the quench at t = 
(completely disordered initial state) and {<7j(i)} is the 
spin configuration at time t. 

First of all we calculate C{t w +t,t w ) in the ferromag- 
netic part of the phase diagram, i.e., for J < Jt. We 
found that C(t w + t,t w ) depends on t and t w through 
the ratio t/t w , as shown in FigfTBl This type of scaling 
is called simple aging and it is characteristic of a simple 
coarsening (i.e., domain growth) process. This result is 
in agreement with the observed behavior of the excess of 
energy (Fig[9|) 

Next we consider the behavior of the correlations dur- 
ing the glassy regime observed in the previous section for 
J > Jt- The typical behavior of the autocorrelation func- 
tion is shown in Figure \17\ where we plot C(t w +t,t w ) vs 
t for T — 0.04, J — 2 and L = 256 and different waiting 
times. The simulation was run up to t — 5 x 10 5 MCS 
and typical averages were performed over 2000 realiza- 
tions of the thermal noise; both times were chosen such 
that ro < t+t w < T\, A first trial to collapse those curves 
showed that in this case C(t w + t,t w ) does not exhibit 
simple aging. A similar behavior is observed for J = 2, 
T = 0.06 < T g and r < t + t w < n. 
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FIG. 17: (Color online) Two-times autocorrelation function 
as a function of the difference of times t for J = 2, T = 0.04, 
L = 256 and different values of the waiting time t w (increasing 
from left to right). 



It has been observed for a large variety of systems that 
in the aging scenario the curves of C(t w + t,t w ) for differ- 
ent t w always collapse into a single one using an adequate 
scaling function^. Although there is no theoretical ba- 
sis for determining the scaling function, there are a few 
choices that have been able to take into account both 
experimental and numerical data, perhaps the most fre- 
quent one being the additive form 

c(t w + 1, t w ) = c st (t) + c ag + _ (10) 



V h(t w ) 



where C s t{t) is a stationary part, usually well described 
by an algebraic decay 



C st (t)=Bt- 



(11) 



The function h, appearing in the aging part of the auto- 
correlations C ag , is some scaling function (in the case of 
simple aging h(t) is a power law that describes the char- 
acteristic linear domain size growth) . In our case the best 
data collapse of the autocorrelation curves was obtained 
using a scaling function of the form 



h(t) = exp 



1-M 



i-m" 



(12) 



which has been used to account for experimental data 23 
and in the Edwards-Anderson model for spin glasses^, 
where r is a microscopic time scale. It is worth to note 
that the scaling function lfl2|) interpolates a range of sce- 
narios: from sub-aging for < \i < 1, to super-aging for 
fi > 1, through simple aging for 23 fi = 1; for /i = TTI is 
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FIG. 18: (Color on-line) Data collapse of the autocorrelation 
curves of Fig[T7](J = 2, T = 0.04 and L = 256) using the scal- 
ing function Eg. (|12[) with scaling parameters shown in Table 
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TABLE I: Scaling parameter values obtained from the best 
data collapse for the correlation curves for J = 2 at different 
temperatures using the scaling forms Eqs. lfT0l) - (fT2|) . 



recovered. In Figure fl8l we show the collapse of the data 
from Figure [17] using the scaling parameter values shown 
in Tabled (r was arbitrarily fixed to one). A similar data 
collapse was observed for J = 2 and T = 0.06 (see param- 
eters in Table [J). To check possible finite size effects we 
performed a similar calculation for T — 0.04, J = 2 and 
L = 64, finding the same collapse shown for L — 256 in 
Figure HU with the same scaling parameters; only a small 
variation in the master curve is observed. We see that 
for T < T g the best data collapse is obtained without 
stationary part and a clear sub-aging is observed. 

We also repeat the correlation calculation for temper- 
atures T > T g at time scales corresponding to the or- 
dered regime t w + t > n(T = 0.6 and T = 1 for J = 2 
and L = 64). Again, aging is observed in this regime 
and a data collapse similar to that shown in Figure fl8l 
using the scalings (fT0ll -lfT2 |) is obtained. The correspond- 
ing scaling parameter values are shown in Table [H We 
see that for this temperature range the best data col- 
lapse is obtained by including a stationary part and that 
the scaling parameter fi decreases systematically as the 
quench temperature increases signaling that the system 
approaches TTI. We find that \i becomes zero at a tem- 
perature T w 1.5 < T c , which can be considered as the 
onset of this non exponential relaxation. 



To further characterize this non equilibrium behavior 
we also analyze the generalized Fluctuation-Dissipation 
Relations (FDR), which can be expressed as^: 



R(t w + t, t w ) — 



X(t w + t,t w )dC(t w +t,t w ) 



T 



dt, 



(13) 



where R(t w + t,t w ) = l/N J2i 9 M*«< + *)) /9hi{t w ) 
is the response to a local external magnetic field hi(t) 
and X(t w +t, t w ) is the fluctuation dissipation factor. In 
equilibrium the Fluctuation Dissipation Theorem (FDT) 
holds and X(t w + t,t w ) — 1, while out of equilibrium 
X depends on t and t w in a non trivial way. It has been 
conjectured 2 ^ that X(t w + t,t w ) = X[C(t w + t, t w )]. This 
conjecture has proved valid in all systems studied to date. 

Instead of considering the response function it is easier 
to analyze the integrated response function 



t m +t 



R{t w + t, s) ds. 



(14) 



Assuming X(t w +t,t w ) — X[C{t w + t, t w )] one obtains 



T X (t w +t,t w )= / X(C)dC (15) 

JC{t w +t,t m ) 

and by plotting Tx vs. C one can extract X from the 
slope of the curve. In particular, if the FDT holds X = 1 
and T\(t) = (1 — C(t)); any departure from this straight 
line brings information about the non-equilibrium pro- 
cess. In numerical simulations of spin glass 27 , structural 
glass2& and Random Anisotropy Heisenberg 2 ^ models it 
has been found that, in the non-equilibrium regime, this 
curve follows another straight line with smaller (in ab- 
solute value) slope when t/t w 3> 1. In this case the 
FD factor X can be interpreted in terms of an effective 
temperature 3 ^ T eff = T/X. 

We apply this procedure during the glassy and ordered 
relaxation regimes previously found for J > Jt. At time 
t w we took a copy of the system spin configuration, to 
which a random magnetic field hi(t) — ha was applied, 
in order to avoid favoring long range orde r 31 i 32 ; a was 
taken from a bimodal distribution (e^ = ±1). We have 
used different values of h in order to check that the sys- 
tem was within the linear response regime. All the results 
presented here were obtained with h = 0.025. 

In Fig. [19] we display Tx(t,t w ) vs C{t w + t,t w ) in a 
parametric plot for different waiting times in the glassy 
regime (i.e., t w +t < n) at T = 0.04. We observe a typ- 
ical two time separation behavior—. At t = the system 
starts in the right bottom corner (fully correlated and 
demagnetized) and during certain time (that depends on 
t w ) it follows the equilibrium straight line, indicating the 
existence of a quasi equilibrium regime. Nevertheless, at 
certain time the system clearly departs from this quasi- 
equilibrium curve and moves along a different straight 
line, but with a different (smaller) slope, indicating an 
effective temperature that is larger than the temperature 




0.7 0.8 
C(t +t,t ) 



FIG. 19: Parametric plot of Tx{t,t w ) vs C(t w +t,t u,) in the 
glassy regime (J = 2, T = 0.04 and L = 256;t w + t < n) for 
different waiting times. The dashed lines are linear fittings; 
the dash-dotted lines represents the thermal equilibrium rela- 
tion T\ = 1 — C. The linear fittings have a slope X = 0.01, 
corresponding to an effective temperature T e // = 4. Waiting 
times from top to bottom are 8 x 10 4 , 10 5 and 2 x 10 5 MCS 
respectively. 



of the thermal bath (T e ff = 4 for the data of FigJTSJ) . 
Notice that the quasi equilibrium regime is very small, 
consistently with the absence of a stationary part ob- 
served in the correlation function (see Table [J) . 

In Fig. [20] we display T\(t,t w ) vs C(t w + t,t w ) in a 
parametric plot for different waiting times in the ordered 
regime (i.e., t w +t > n) at T — 0.6. The two-slope struc- 
ture is again observed, although with a smaller effective 
temperature (T e // = 1.5). 



IV. DISCUSSION 

We introduced a lattice spin model that mimics a sys- 
tem of interacting particles through a short range repul- 
sive potential and a long range power law decaying poten- 
tial. Through a detailed Monte Carlo simulation analysis 
we computed the complete equilibrium phase diagram of 
the model at finite temperature and characterized the or- 
der of the different transition lines. We showed that the 
model presents only two simple ordered phases at low 
temperatures: ferromagnetic (for J < J<) and antiferro- 
magnetic (for J > J t ), without any trace of geometrical 
frustration and/or complex patterns. 

We then analyzed the out of equilibrium relaxation 
of the system after a quench from infinite temperature 
down to subcritical temperatures, in different regions of 
the phase diagram. While a normal coarsening behavior 
appeared in the ferromagnetic region of the phase dia- 
gram J < Jt (i.e., a domain growth process that follows 
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FIG. 20: Parametric plot of Tx(t,t w ) vs C(t w +t,t w ) in the 
ordered regime (J = 2, T = 0.6 and L = 64; t w + 1 > n) for 
different waiting times. The dashed lines are linear fittings; 
the dash-dotted lines represents the thermal equilibrium re- 
lation Tx = 1 — C. The linear fittings have a slope X = 0.4, 
corresponding to an effective temperature T e // = 1.5. Wait- 
ing times from top to bottom are 4 x 10 4 , 8 x 10 4 and 10 5 
MCS respectively. 



Allen-Cahn law lit) ~ t x l 2 ) the system shows a com- 
plex relaxation scenario in the antiferromagnetic region 
J > Jt- This is precisely the most interesting situation, 
since for those values of J the pair interaction potential 
of the model shows the same qualitative features as a 
continuous, Lennard-Jones (LJ) like potential, namely, 
a nearest neighbors repulsive interaction (i.e, "hard core" 
like) and an attractive, power law decaying interaction at 
longer distances (see FigfTJ) . We must stress that we did 
not intend to present a lattice version of the LJ gas, but 
to show that some very basic features present in it (i.e., 
the competition between short and long range interac- 
tions) are enough to produce non trivial slow relaxation 
properties. We observed that such competition gives rise 
to non-equilibrium structures that strongly slow down 
the dynamics, even in the absence of geometrical frus- 
tration and/or imposed disorder. The most interesting 
of those structures give rise to a relaxation regime with 
several appealing properties that strongly resemble those 
observed in different glassy systems. 

First of all, that regime is characterized by a dynam- 
ically generated disordered non-equilibrium state, char- 
acterized by a labyrinth structure, i.e., composed of at 
least one macroscopic connected domain. The energy of 
such state displays a pseudo-plateau and a finite lifetime 
Ti that diverges when T — > 0; both properties are in- 
dependent of the system size. Such phenomenology is 
extremely reminiscent of transient particles colloidal gels 
obtained by quenching a monodisperse gas of colloidal 
particles under Brownian Dynamics (Langevin dynamics 



in the overdamped limit) that interact through a general- 
ized 2n — n LJ potential o 34 ' 35 . A gel is a non-equilibrium 
disordered state, characterized as a percolating cluster of 
dense regions of particles with voids that coarsen up to 
certain size and freeze when the gel is formed; transient 
gels do not have permanent bonds between them and col- 
lapse after a finite life time 3 ^. The energy of transient gels 
of 2n — n LJ colloidal particles displays a slowly decaying 
pseudo plateau 3 - 6 -, as in the present case. It has also been 
observed that the life time of the gel strongly increases 
as the interaction range is decreased 36 , for instance, by 
changing the value of n. It would be interesting to check 
if a similar effect can be obtained in the present model by 
changing the interaction range (for instance, by chang- 
ing the exponent of the long range interaction term in 
Hamiltonian ([2])). 

Second, during the glassy regime the dynamics appears 
to be governed by free energy barriers to coarsening that 
scale as a power law with the characteristic domain size, 
with an associated logarithmic growth l(t) ~ [T ln(t)] a 
(and therefore a logarithmic relaxation). Such behavior 
is characteristic of class 4 systems, according to Lai et al 
classification? . Some examples of non-disordered short- 
range interacting systems that present growing free en- 
ergy barriers with the domain size are already known, 
such as the three dimensional Shore and Sethna (SS) 
mode l 18 i 37 i 38 (i.e., an Ising model with nearest neigh- 
bors ferromagnetic interactions and next nearest neigh- 
bors antiferromagnetic interactions) and a generalization 
of the previous one introduced by Lipowski et in- 
cluding a four spin plaquette interaction term. However, 
in those models the barriers appears to grow linearly with 
the domain size (which corresponds to a pure logarith- 
mic growth l(t) ~ ln(i)) and therefore fall into the class 
3 category of Lai et aU£. So far, examples of class 4 sys- 
tems were found only among disordered systems, such as 
spin glasses^ and the Ising model with random quenched 
impurities^ -. To the best of our knowledge, this would 
be the first possible example of class 4 behavior in a non 
disordered system. Nevertheless, an important difference 
between the above mentioned non-disordered models and 
the present one have to be remarked. While in those 
models the logarithmic growth appears to lead to diver- 
gent barriers, in the present case the apparently barriers 
growth stops at some maximum value F that determines 
the life time t\ . This implies the existence of a character- 
istic length Imax in the dynamics of the system. Although 
such limited length scales makes very difficult to obtain 
better numerical evidence of the existence of growing bar- 
riers, the consistency between the scaling of the excess of 
energy and the time for shrinking square excitations gives 
support to our conjecture. Thus this system appears to 
behave, at least for certain time scale that can become 
very long a very low temperatures, as a class 4 system, 
even though in the long term it behaves as a class 2 sys- 
tem in the sense that its ultimate dynamics is governed 
by a single free energy barrier. This opens the possibil- 
ity of having a truly non disordered class 4 system if the 
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lifetime of the glassy state at finite temperature could 
be extended by tuning the range of the interactions, as 
previously discussed. Indeed, we believe that the possi- 
bility of having in a non-disordered class 4 system makes 
it worth to further investigation. 

While we do not have an explanation for such possi- 
ble relaxation scenario (i.e., power law growth of barriers 
with a rather small exponent l/a and limited length scale 
for growing) probably a key ingredient to explain it would 
be the moderated long range character of the interac- 
tions. It would be very interesting to check if the same be- 
havior can be detected in the LJ gas or its generalizations, 
which appear to present a similar phenomenology 36 . 
However, the range interactions would be not enough in 
the present model to generate dynamical frustration (in 
the sense stated by Shore et at!£, i.e., systems whose dy- 
namics is slowed down by the presence of growing free 
energy barriers), but the type of competition between 
interaction would be equally important. This can be 
clearly seen by looking at the non-equilibrium behavior 
after a quench of the reverse model, namely, that given 
by the Hamiltonian (U) with the inverse coefficients sign 
(Ji < and J2 > 0). While the equilibrium phase di- 
agram of that model is by far more complex than the 
present one&£, its domain growth behavior after a quench 
to subcritical temperatures is relatively simple, at least 
for the regions of the phase diagram explored up to now. 
Depending on the ratio of couplings J\ j J2 it behaves as 
a class 2 system^ (which implies domain-size indepen- 
dent free energy barriers), like the 2D SS mode l 18 ' 37 , 
or relaxation can be dominated by nucleation effects 4 ^. 



Accordingly, that system presents simple agin g 43 ' 44 and 
trivial FD relations 3 -^ (infinite effective temperature), at 
variance with the present case. 

In the glassy regime of the present model we found non 
trivial aging effects, with scaling properties characteris- 
tic of glassy systems (subaging). Non trivial aging has 
also been found in the non disordered four-spin ferromag- 
netic model 4 ^ (a particular case of the model of Lipowski 
et aU^), but in this case the system displays superaging, 
while the disordered version of the same model displays 
subaging 4 ^ Subaging has also been reported in molecu- 
lar dynamic simulations of small LJ clusters^ 7 -. 

We also found FD relations displaying a well defined ef- 
fective temperature in the aging regime. It is interesting 
to note that, although non trivial behavior of time cor- 
relations and responses are usually associated to glassy 
systems, they have also been observed in very simple sys- 
tems which undergo domain growth at intermediate time 
scales (as in the present case) , namely, the ferromagnetic 
Ising chain 4 ^ and the 2D ferromagnetic Ising and Potts 
models with Kawasaki dynamics^. A similar behavior 
has been found in a non-disordered plaquette model for 
glasses introduced by Cavagna et a l 49 i 50 . It is worth not- 
ing that the 3D SS model presents only trivial FD rela- 
tions, even for the temperature range where it appears 
to present logarithmic growth of domains 3 ^. 
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